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radiation is beamed into a wide but hollow conical radiation pattern. The 
emission then reflects the rotation of the source in that a remote observer 
first detects radiation from one side of the emission cone and then the other 
[for a review, see Smith, 1976; Carr et al., 1982; Goldstein and Goertz, 
1982], Dulk [1965] concluded that an emission cone angle of 79° worked well 
for the lo controlled sources. Goldstein et al. [19793 came to similar 
conclusions for the lo independent sources. Because radiation emitted from 
the back side of Jupiter cannot reach Earth, one would expect that the phase 
of lo would always be greater than 90°, but less than 270*. In fact, radia- 
tion from Io-6 is often received at lo phases less than 90°. This introduces 
an unexplained asymmetry in the occurrence probability of the lo phase depend- 
ent source. In this paper we propose a new explanation of this asymmetry. 

The asymmetry has been known for some time, and several authors have 
offered a variety of explanations for it. Dossier and Hill [19793, following 
earlier ideas [see, e. g., Goldreich and Lynden-Bell, 19693, pointed out that 
the difference between the sub-Io longitude and the longitude of the footpoint 
of the lo flux tube (IFT) is about 30° when the sub-Io longitude is 230° 
(using the GSFC/Ojj model for the global Jovian magnetic field [Acuna and Ness, 
19763), which is about what is needed. However this effect is in fact can- 
celed by the rather large azimuthal deviations of the magnetic field from the 
meridional plane in the relevant sub-Io longitudes. Piddington and Drake 
[19683 considered the bending of the magnetic field by lo and estimated the 
propagation time of hydrcmagnetic waves from lo to the Jovian ionosphere. 
They assumed a velocity of propagation along the lo flux tube of 500 km/s, or 
a propagation time of about 20 minutes. According to recent models a^'duoed 
from Pioneer 10 and 11 and Voyager 1 and 2 data, the Alfven velocity along the 
IFT is almost equal to the velocity of light, except in the lo torus and in 
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ABSTRACT 


We propose sn explanation of an asymmetry in the occurrence 
probability of the lo-^ependent Jovian decanetric radiation, 
lo generates stronger Alfven waves toward the south irtien it 
is in the northern part of the torus. This wave then gener- 
ates decametrlc radiation in the northern ionosphere after it 
reflects in the southern ionosphere. The asyraaetry then re- 
sults from computing the propagation time of the Alfven wave 
along this trajectory. The ray paths of the decameter radia- 
tion are calculated usir% a three-dimensional ray tracing 
program in the Jovian ionosphere. Variations in the expected 
probability plots are computed for two models of the Jovian 
ionosphere and global magnetic field, as well as for several 
choices of the ratio of the radiated frequency to the X-mode 
cutoff frequency. 
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1. Introduction 


The Jovlen deoametric radiation, or DAH, extenda in frequency fron below a 
few negahertx to nearly 40 HHz. It la well known that the probability that a 
remote obaerver will receive radiation from Jupiter at decameter wavelengtha 
ia a critical function of the central meridian longitude (CML) of the planet 
aa aeen by the obaerver. In 1964 Bigg discovered that the phase of the moon, 
lo, controls the occurrence probability of the most intense decameter storms 
do phase is defined as the angular departure of lo from superior conjunction 
as seen by an observer). A second component is independent of lo phase, but 
ia observed at central meridian longitudes close to those of the lo phase 
dependent component. The occurrence probability of DAM is often shown in two 
dimensional plots in which the axes are central meridian longitude (in System 
III) and lo phase. The occurrence probability is then indicated by varying 
shades of grey or by contours of increasing probability. The most prominent 
features in the occurrence probability plots are the lo phase dependent 
sources; the main source (lo-A, phase of lo ~ 240* and CML - 230*) and the 
early source (lo-B, phase of lo - 90* and CML - 150"). In addition, there are 
lo phase independent features at nearly the same CML’s as for the lo A and B 
sources as well as an additional lo-dependent source at phase of lo - 240" and 
CML - 300“ (source C). Other minor features exist v^ich are described in some 
detail in the recent extensive review of decameter phenomenology by Carr et 
al. [1982]. 

One of the most appealing explanations for the observation that the deca- 
meter radiation comes primarily from two regions in CML and from two positions 
in lo's orbital phase is the idea first proposed by Dulk [1965] that decameter 


radiation is beamed into a wide but hollow conical radiation pattern. The 
emission then reflects the rotation of the source in that a remote observer 
first detects radiation from one side of the emission cone and then the other 
[for a review* see Smith. 1976; Carr et al., 1982; Goldstein and Goertz, 
1982]. Dulk [1965] concluded that an emission cone angle of 79** worked well 
for the lo controlled sources. Goldstein et al. [1979] came to similar 
conclusions for the lo independent sources. Because radiation emitted from 
the back side of Jupiter cannot reach Earth, one would expect that tte phase 
of lo would always be greater than 90°, but less than 270°. In fact, radia- 
tion from lo-B is often received at lo phases less than 90°. This introduces 
an unexplained asymmetry in the occurrence probability of the lo phase depend- 
ent source. In this paper we propose a new explanation of this asymmetry. 

The asymmetry has been known for some time, and several authors have 
offered a variety of explanations for it. Dessler and Hill [1979], following 
earlier ideas [see, e. g., Goldreich and Lynden-Bell, 19693, pointed out that 
the difference between the sub-Io longitude and the longitude of the footpoint 
of the lo flux tube (IFT) is about 30° when the sub-Io longitude is 230° 
(using the GSFC/Ojj model for the global Jovian magnetic field [Acuna and Ness, 
1976]), which is about what is needed. However this effect is in fact can- 
celed by the rather large azimuthal deviations of the magnetic field from the 
meridional plane in the relevant sub-Io longitudes. Piddington and Drake 
[1968] considered the bending of the magnetic field by lo and estimated the 
propagation time of hydrcmagnetic waves from lo to the Jovian ionosphere. 
They assumed a velocity of propagation along the lo flux tube of 500 km/s, or 
a propagation time of about 20 minutes. According to recent models aoduoed 
from Pioneer 10 and 11 and Voyager 1 and 2 data, the Alfven velocity aiot^g the 
IFT is almost equal to the velocity of light, except in the lo torus and in 
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the Jovian ionosphere. Furthermore, when lo is in the northern-nost part of 
the torus, the effect of the bending of the aagMtio field lines is almost 
negligible. McCulloch [1971] performed ray tracing of irtiistler mode raves 
from lo to the ionosphere and calculated both propagation times and the 
relative intensities to arrival points in the "coupling region" using 
Melrose's [19671 magnetospheric density model. McCulloch did not consider any 
generation mechanisms for the whistler mode waves. Furthermore, the density 
model that he used is not supported by recent observations. Goertz and Delft 
[1973] assumed that an Alfven wave reflected in the ionosphere reaches lo in a 
short time and doubles the magnetic tilt. They suggested that the asymmetry 
could arise from reconnection at an X-type neutral point in a plasma wake of 
lo. Smith [1973] assumed that the plasma instability that produces DAM could 
be convective so that in the source region very small values of the ip'oup 
velocity could lead to long delays (more than 10 minutes would be required) in 
the escape of the radiation into the magnetosphere. 

The most detailed comparison between the observed occurrence probability 
and the assumed beaming pattern of the lo dependent source has been carried 
out by Thieman and Smith [1979]. They assumed that radiation emerged in a 
thin conical sheet centered on the magnetic field line at the point of its 
intersection with the Jovian ionosphere. Any reflection or refraction caused 
by the ambient plasma was neglected, but they did include occultatlon by the 
planet. They found that if the source region was located alor« a field line 
20" ahead of lo (the "lead angle") the asymmetry could be accounted for. The 
best agreement with observations was obtained using a cone angle of 80*. In 
this paper we propose an explanation of this lead angle in which we Include 
the effects of wave refraction in the Jovian ionosphere and incorporate a 
model of Alfven wave propagation through the lo plasma torus. The basis for 


this model is given in section 2. Three dimensional ray tracing is performed 
in order to explain the occurrence probability map and to evaluate conditions 
for wave generation in the Jovian ionosphere including constraints on the wave 
normal angles and on the ratio of the generated frequency to the right hand 
cutoff frequency. The ray tracing calculation is explained in section 3 and 
the results of the calculation are compared with observations in section 
In section 5, we present our conclusions, 

2. lo phase asymmetry 

Neubauer [1980] has shown that because of the existence of the lo plasma 
torus, the unipolar current between lo and the Jovian ionosphere described by 
Goldreich and Lynden-Bell [19693 is in fact carried by large amplitude Alfven 
waves (or Alfven current wings). It has been suggested by many authors 

[Gurnett and Goertz, 1981; Goldstein et al., 1982; among others] that the 

strong parallel electric fields in these "kinetic” Alfven waves [Hasegawa and 
Chen, 1976; Hasegawa and Mima, 1978] are the origin of the keV electrons that 
are thought to generate DAM. 

One of the striking observations made by the Planetary Radio Astronomy 

(PRA) experiment on Voyagers 1 and 2 was the discovery that most of the 
decameter radiation has a dynamic spectrum which resembles nested parentheses, 
or arcs, on a frequency-time plot [Warwick et al., 1979]. The arc pattern was 
shown to be consistent :vith Dulk's [1965] emission cone model by Pearce 
[1981], Staelin [1981], and Goldstein and Thieman [1981]. Gurnett and Goertz 
[1981] proposed that the discreteness of the arc pattern could be a result of 
multiple reflections of the standing Alfven wave system [but also see 
Lecacheux et al., 1981]. The propagation time of the Alfven waves has 
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recently been calculated by Bagenal [1982], and we make detailed use of her 
calculations in the discussion below. The intensity of the Alfven waves 
generated by a satellite was calculated in a linear treatment by Drell et al. 
[ 1965 ]. Neubauer's [1980] work is its extension to the nonlinear (large 
amplitude) regime [for reviews, see Goldstein and Goertz, 1982; Hill et al., 
1982]. It was found that both the magnetic field intensity and the power in 
the generated Alfv8n waves are approximately inversely proportional to the 
Alfven velocity, c B^//(Mirp) where p is tl» ion mass density and is the 

magnetic field [see Neubauer, 1980 for details], 1. e., proportional to the 
square root of the mass density. 

These calculations assumed that the density profile surrounding the 
satellite was constant. The wavelength of the waves perpendicular to the 

ambient magnetic field is almost equal to the diameter of lo (2 R]^q, where 

Rj^ is the radius of lo). The parallel wavelength then becomes 

where is the Alfven Mach number. Results from the Voyager 1 Plasma Science 
Experiment indicate that > 0.25 near the peak density of the torus 
[Bagenal, 1982]. Therefore, 

X„ - 8 Rj^ - Rj/5 (2) 

where Rj is the radius of Jupiter. In less dense regions, M^ is smaller and 
X„ is larger. Because the radius of the torus is almost equal to that of 
Jupiter, the density profile surrounding lo cannot be assumed constant. The 
torus density surrounding lo is significantly smaller than its peak value when 
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lo reaches Ita most distant excursions north and south of the Jovian nagmtic 
equator [see. Warwick et al., 1979!!> 17te Alfven wave intensity is controlled 
by the mean density over at least several wavelengths. When lo is not in the 
center of the torus, the mean densities over several wavelengths differ consid- 
erably in the northern and southern parts of lo's excursion in the plasma 
torus because the wavelengths tend to be larger in the lower density regions. 
Because of this, the Alfven wave intensities should differ to the north and 
south of lo. When lo is in the northern-most part of the torus, where strong 

DAH is usually received, lo sees a denser plasma in the south and a dilute 

plasma in the north. In this situation a stronger wave is generated south- 
ward. Because the northward directed wave is weak, it is unlikely to stimu- 
late decameter radiation when it reaches the northern ionosphere. However, 
the intense southward directed wave will, after reflection in the southern 
ionosphere, be able to generate DAM in the northern ionosphere. Gurnett and 
Goertz [1981] estimated that absorption of these Alfven waves was sufficietitly 
weak that tens of reflections were possible before significant attenuation 
occurred. Let us leave aside the question of whether observable radiation 
will be generated in the southern hemisphere and assume that all the observed 
radiation originates in the northern ionosphere. This assumption is 
consistent with the PRA observations that most of the decameter radiation 

appears to originate in the northern hemisphere [see Carr et al., 1982], In 

that case it is possible to evaluate the expected asyranetry if radiation 
occurs only after a reflection in the southern ionosphere. 

The propagation time of an Alfven wave using the WKB approximation has 
been calculated by Bagenal [’982] (see Fig. 7 of that paper). Her results, 
which may require some modification due to large density gradients in the 
torus, can be approximated as follows; 
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T do to southern ionosphere) > 6.5 3.5 sin # (min.) 

9 

<3) 

do to northern ionosphere) ■ 7.0 - 4.3 sin * (min.) 

where • > e^^Csub-Io longitude) 112*. The propagation time for the wave to 
reach the northern ionosphere after reflection in the southern ionosphere T is 

T s 2T. ♦ a 20 ♦ 2.7 sin ♦ (rain.) (4) 

9 n 

lo moves 0. 14* /min and Jupiter rotates 0.6*/min, so that the apparent phase of 
lo changes 0. 46T* during the propagation time. The delay angle during T 

minutes then becomes 

a 9.6* ♦ 1.24* sin « (5) 

During its two traverses of the torus the wave is significantly refracted by 
an angle 

= 40’M^ (1.5 ♦ 0.5 sin •) (6) 

where is effectively 0.125 according to Fig. 8 of Bagenal [1982]. Then 

Ae^ r 7.5* ♦ 2.5* sin • (7) 

The total lead angle is an addition of (5) and (7). or 
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« 17.1* ♦ 3.7«* sin • 


( 8 ) 


This value is close to the 20* found by Thieman and Smith [19793 needed to 
give the best agreement with the observed occurrence probabilities. We now 
turn to consider the role of refraction of the decametric X-roode radiation in 
the ionosphere. 


3. Three*>dimenslonal Ray Tracing 

It is commonly, though not universally, assumed that DAM is generated just 
above the R mode cutoff frequency f^ (i. e., the cutoff frequency of the right 
circularly polarized fast extraordinary mode, see, e. g., Stix [1962]) at 
large wave normal angles and then propagates to Earth in the extraordinary (X) 
mode [see the recent review by Goldstein and Goertz, 1982], The frequencies 
must be higher than fp in order that DAM can escape from Jupiter in the R-X 
mode. When the frequency of these waves is close to fp, they will be re- 
fracted toward the direction of the magnetic field and thus the effective cone 
angle of the emission pattern becomes narrower. Waves propagating toward the 
planet are reflected around the ionosphere. In order to evaluate such plasma 
effects, three-dimensional ray tracing was performed for waves generated at 
various locations around the ionosphere. The numerical program is based on 
one prepared for ray tracing in the ionophere of Earth by Jones and Stephenson 
[1975], and modified for use with Jovian parameters. 

In using this program it is necessary to specify a global model of the 
Jovian magnetic field. Most of our calculations were done using the GSFC/0^ 
model [Acuna and Ness, 1976; Acuna et al., 1982], but where difference" arc 
significant, the 16 coefficient JPL octopole model with three internal terms 
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and on# axtarnal tarn (13^1) was usad for oonpariaon CSnith at al., 1976]. 
Both of thaaa nodala ara basad on Pionaar 11 data. In tha aagnatosphara » tha 
alaotron dansity distribution was astinatad using Santaan and Goartz's [1978] 
alaotron dansity aodal. Any af facts of ions wara naglaotad. Tha ray tracing 
calculation included the oblatanass of Jupitar. Tha radial distance to Um 
cloud tops R (in ka) war taken to be 

R * 71372//(1 - 0.14371 cos^e) (9) 

where « is colatitude. We have slightly modified Sentaan and (!oartz*8 density 
model In that ye use R as defined in (9) in place of their Rj. Whan tha radia- 
tion is propagatif^ in the topside ionosphere, as it does near tha generation 
region, than a model for the ionospheric electron distribution must be usad. 
Atraya and Donahue [1976] have derived such a modal basad on Pionaar 10 obser- 
vations. We have used a nonlinear least squares fit to the results of their 
modal in order to obtain an analytical approximation of their results. Tha 
approximation wa use is 

4 ^1 

5 h - 1417 991 
"l “ ’••'5 » 366.3 > * <r^> = 

where N^(cm ) is the ionospheric electron density, h is the altitude in km 
above tha cloud tops and s 0 for h < 276 km. By combining the 
contributions to the ionosphere model with that of the magnetosphara model wa 
obtained a single electron density modal to use in the refraction calcula- 
tions. This combination is included in the fitting. The model is shewn as 
the solid curve in Figure 1. Tte effects of the plasma torus are not included 
because the dansity is too low to affect the ray paths of DAM. because the 
maximun plasma frequency in the torus is only 0.5 MHz [Warwick at al.. 1979J. 
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The waves are launched under the following starting conditions: 


1. The frequency f is set to some fraction above the local R mode cutoff 
frequency fj^ along the lo flux tube. This fraction is varied from 0.19 to 

109. 

2 . Initial wave normal angles of 70*, 80*, and 90* are used. 

3. The conical radiation pattern is simulated by initializing the azimuthal 
angle in 10* steps about the magnetic field. 

Ray paths are calculated up to an altitude of 2.5 Jovian radii above the 
cloud tops where it is then assumed that plasma effects are negligible. The 
final ray direction is defined by two angles: The angle of the r?y with respect 
to the zenographic equatorial plane (d) and the equivalent longitudinal 
direction A, the angle between the projection onto the equatorial plane of the 
ray direction and the direction of CML = 0*. Rays with |d| <2* are assumed 
to reach Earth. DAM is received when the line of sight to the observer is 
along A. 


4, Comparison between ray tracing results and observations 

The following calculations were performed under the assumption that DAM 
occurs only In the northern hemisphere and then only when lo is in the 
northern part of the torus. 

Figure 2 is an example of a calculation for f s 22 HHz and 1.01 (19 

above fp). Results are shown only for sub-Io longitudes between 120* and 270“ 
where lo is in the northern part of the torus; 30“ steps in sub-Io longitude 
were used and ray directions are computed only for a northern hemisphere 
source. The dashed line Indicates 120* sub-Io longitude. The symbols X, Y, 
and Z correspond to rays Initially launched with wave normal angles ol 7Cr , 
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80* t and 90*. raapaotively. The ayart>ola are indioated only if |d| < while 
the ray traoli^ wea perforned with 36 initial aiimtthal anglea. When )d) is 
between 2* and the sizes of the synbols are chained continuously (larger 

syabols correspond to saalJer values of |d|). When |d| <2* the line widths 

of the syabols are systeBatically varied in four discrete steps so that 

thicker symbols are used for smaller values of |d|. The phase of lo (yj^) is 

calculated by the equation: 

Yio s ISO* CML - - &♦ (11) 

where Is sub-lo longitude and &• is the lead angle as defined by (8). 

Notice that there is no contribution from sub-Io longitude 270* because there 
the initial condition that f * 1.01 fj^ cannot be satisfied; the required iono- 
spheric density is too low at the altitude corresponding to a 22 MHz gyro- 
frequency. In fact from (10), because the required altitude is h < 276 km, 
the density given by the model is zero. The solid curves in the figure show 
contoia* plots of occurrence probability at 22 MHz based on Thieman (1979). 
The contours are slightly different from Thieman and Smith's [1979) because of 
a different smoothing method. In this figure a® * 0*, that is, no lead angle 
is assumed. This model, of course, includes the difference between sub-Io 
longitude and the longitude of the IFT foot point as well as the declination 
of the magnetic field using the 0^ field. Note that, as expected, lo phases 
less than 90* are not allowed, and the results do not fit the observations. 

Figure 3 includes the lead angle calculated from (8). Other parameters 
are the same as those of Figure 2. I*"' phases less than 90* are now apparent. 

Fairly good agreement with observations is seen especially for wave normal 
angles 60* and 90*. Radiation is also allowed from CML ~ 300*, lo-C. Note 
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that the lower part of the lo-A contour here and in Figure 2 includes radia- 
tion ft*om the lo-lndependent source region and is not applicable to this 
model. 

In Figure 4 we show the result of assuming that f = 1.001 f^ (0.1% above 

f|^) at 22 MHz. Now the radiation frequency is so close to f^ that increased 

refraction causes the effective cone angle to become smaller. Any generation 
mechanism which predicts that the radiation frequency is this close to 
cannot account for the occurrence probability of DAM so long as one retains 
the basic hypothesis that the A and B source regions represent emission from 
two sides of the same conical radiation pattern. At higher ratios of f/f|^. 
for example f = 1. 1 fj^, (10% above fj^), even better agreement with observa- 
tions (cf. Figure 3) is found, especially when fairly small wave normal angles 

(~ 70*) are used (Figure 5). However, the results corresponding to sub-Io lon- 
gitude near 270" do not fit wll. The source region in this longitude range 

is at a very low altitude, just above the peak in the ionospheric density. 

Furthermore, the southern edge of the radiation cone is directed toward the ob- 
server, and thus the angle Id I is very sensitive to the detailed behavior of 
both the magnetic field and electron density models. This choice of f/f^ 

gives an especially good fit in the lo-C source region. This then provides a 
possible explanation for the predominantly right circular polarization ob- 
served from this source in the PRA data [e.g., Carr et al., 1982]. 

The role of refraction is emphasized in Figure 6 , where, for comparison 
with the earlier results (Figure 3), «e have artificially set the radiating 

frequency to 99 MHz, far above fp. Refraction by the plasma is neglected and 

ray paths are frequency-independent in this case. The remainder of the param- 
eters are identical to the case considered in p’igure J. In this case, rays 
can be reflected at the cloud tops, but such rays do not reach the obse‘~vcr . 
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In this way the oooultation Included in ThleMn and Siaith*a Cl 9793 caleula- 
tiona can be reproduced. This figure is alnost identical to their Pig<s*e S. 
In this case, the observations are better fit for aaallw* wave noraal ai^llas. 

We have tried to investigate the sensitivity of our results to our 
particular choice of ionospheric density nodel by repeating sons of the 
calculations using a model developed from Voyager 1 data by Strobal and Atreya 
[1962]. To include the generally higher densities found in this models we 

have increased by a factor of three the magnetospherio density derived by 
Sentman and Goerts [1978]. This model is also shown in the dashed curve of 
Figure 1. The results are shown in Figure 7 for t * 1.01 f^. The higher 
density causes more refraction which reduces the effective cone angle. In the 
case of f s 1.1 f|j, the result is almost same as was shown in Figure 5. 

For the remainder of this discussion we return to using the Pioneer 10 

density model and examine results for other frequencies. In Figure 8 we show 
a calculation using f = 1.01 f^ at 15 MHz. Very good agreement with the ob> 
servations is seen. In this frequency range, much of the radiation from CNL's 
greater than 220* is dominated by lo phase independent emission and cannot be 
treated here. 

In Figure 9, we repeat the calculation at 30 MHz. The agreement with the 
observed occurrence probability is not so good. However, in a case of f s 1.1 
f|^ (Figure 10), fairly good agreement in lo-B source is obtained for wave 
normal angles of 70* and 80*. Discrepancy in lo-^ source may be due to genera- 
tion mechanism. At such high frequencies (and low altitudes) the results are 

very dependent on high order raultipole contributions to the magnetic field. 
Differences between the various global field nwxlels become noticeable. In 
Figure 11 we have repeated the calculation at 30 MHz using the JPL octopole 
model for f « 1.01 f|^. The agreement with the occurrence probability data is 
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better in some respects in this case than that for f « 1. 1 f|^. At still 
higher frequencies (~ 38 MHz) none of the magnetic field models appear capable 
of correctly locating the lo^ source. This has been a problem with all 
models of the field and has been discussed in the review by Goldstein and 
Goertz [1982]. 

5. Conclusions 

We have proposed a new model to explain the observed lo phase asymmetry of 
the lo controlled A, B, and C decameter sources. In this model we have 
assumed that when lo is in the northern part of the torus an intense Alfven 
wave is generated which propagates southward, reflects from the southern 
ionosphere, recrosses the torus and then generates DAM when it encounters the 
northern ionosphere. The other Alfven wave, initially propagating northward, 
encounters a lower average electron density in the torus which reduces its 
intensity so that it will be unlikely to generate DAM when it reaches the 
northern ionosphere. Standing Alfven waves were observed by the Voyager 1 
magnetic field experiment in the southern hemisphere when lo was in the south- 
ern part of the torus (CML - 330“ ) [Acuria et al,, 19813. According to our 
theory, stronger waves can be expected in the northern hemisphere at that 
time. Our model suggests that the Alfven wing which stimulates the DAM will 
have a higher current density than that deduced by AcuRa et al. [1981]. 
However, when Voyager was close to Jupiter, the PRA experiment did detect DAM 
at nearly all central meridian longitudes [Carr et al., 19823, which indicates 
that the effect we arc describing is probably one of degree and one which 
determines the intensity of DAM rather than whether any radiation is generated 
at all. The fact that DAM seems to be produced predominantly when lo is in 
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the northern half of the torus, 1. e., when the tip of the nurtbern Mnaetiu 
dipole axis points toward lo, is thought to reflect the details of the eals- 
sion process and goes beyond the scope of this discussion. For the sane 
reason, in the context of this nodel, we cannot explain tdiy rMliation is less 
likely to be observed when lo is near the center of the torus, nor can this 
model account for the difference between the maximum radiation fre<|uenoies of 
the lo A and B sources. A related question for which we have no answer, is 
why, when lo is in the southern part of the torus and launotes an intense 
Alfven wave northward (or southward }. DAN is not likely to be observed. These 
problems are probably intimately related to the plasma emission mechanism 
rather than the wave propagation. Finally, the location of the lo^ source at 
high frequencies is not particularly well explained in this model. 

This paper also presents the first three-dimensional ray tracing calcula- 
tions of decameter radiation in the Jovian ionosphere. As a result, we are 
better able to provide realistic limits on some of the source conditions neces- 
sary for the generation of DAM. The best agreement between the ray tracing 
and the observed occurrence probability plots was found for initial wave 
normal angles between 80“ and 90“ and radiation frequencies one to ten per 
cent above fj^. As f/fp increases, smaller initial wave normal angles give 
better fits. The upper limit on f/fj^ cannot be determined from this analysis, 
but will have to be calculated frcm whatever generation mechanism proves 
responsible for (»*oducing DAM. Values of f/fp i 1.001 could not adequately 
model the observations. The density and magnetic field models used in this 
study imply that in the decameter source region the ratio of the plasma fre- 
quency to the cyclotron frequency ranges between 0.02 and 0.09, which is compa- 
rable to some estimates of this ratio in the terrestrial Auroral Kilometric 
Radiation (AKR) source [Calvert, 1981]. These results are consistent with 
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recent theories on generation of DAM and AKR [Wu et al., 1982; Goldstein et 
al., 1982; and the review by Goldstein and Goertz, 1982]. Besides the prob- 
lems mentioned above, further work needs to be done in developing a detailed 
theory of the generation and saturation of the Alfven wings in the presence of 
an inhomogeneous medium. Use of the WKB approximation will probably not be 
appropriate in such an analysis. 
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Figure Captions 

Fig. 1. Jovian electron density models. Squares and circles indicate points 
adopted from the models based on Pioneer 10 [Atreya and Donahue, 1976] and 
Voyager 1 [Strobel and Atreya, 1982], respectively. The solid curve shows 
a Pioneer 10 density model used in our ray tracing based on the data shown 
as squares. The dashed curve shows a Voyager 1 model cased on ^e data 
plotted as circles. 

Fig. 2. The calculated reception pattern in CHL>lo phase plane for 22 MHz and 
f s 1.01 fj^. Initial wave normal angles are 70" (X), 80" (Y), and 90" 
(Z). The GSFC/Oj^ magnetic field model and the Pioneer 10 electron density 
model are used. Contour plots shows smoothed occurrence probability for 
the same frequency based on Thieman [1979]. The contour interval for this 
and all subsequent contour plots is 0.2 starting at 0.1. 

Fig. 3. Identical with Figure 2, but the lead angle defined by (8) is assumed 
in this and the following figures. 

Fig. 4. Identical to Figure 3 except that f = 1.001 fj^. 

Fig. 5. Identical to Figure 3 except that f s 1.1 fp. 

Fig. 6. Identical to Figure 3, but refraction is neglected. Ray paths are 
independent of frequency. Contours are the probablity tar t s 22 MHz. 

Fig. 7. Identical to Figure 3, but the Voyager 1 density model is used. 

Fig. 8. Identical to Figure 3 except that f s 15 MHz. 

Fig. 9. Identical to Figure 3 except that f = 30 MHz (f * 1.01 fj^). 

Fig. 10. Identical to Figure 9 except that f s 1.1 fp. 

Fig. 11. Identical to Figure 9, but the JPL magnetic field model is used. 
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